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ON THE IMPROBABILITY OF FINDING ISOLATED SHOALS IN 
THE OPEN SEA BY SAILING OVER THE GEOGRAPHICAL 
POSITIONS IN WHICH THEY ARE CHARTED. 


By Mr. G. W. Lirrienaes, Washington, D. C. 


Many of the isolated shoals that are represented on nautical charts of the 
oceans have been located from the reports of mariners who have discovered 
them incidentally in making voyages of commerce. Previous to the year 1860, 
when there was no exact knowledge of the depths of the oceans, the vague 
reports of navigators, often doubtless based upon the observation of floating 
objects and of misleading appearances of the surface of the sea, caused the 
charting of many dangers for the existence of which there is no substantial 
foundation. But, as our knowledge of bathymetry increased, the existence of 
many of them was disproved, and they were removed from the charts. 

As a result of these experiences, there arose a traditional distrust among 
mariners and hydrographers of the existence of many of these dangers that 
still appear on the charts with well founded evidence, and there is perhaps a 
disposition on the part of many to claim that they should be removed upon 
scant evidence of their non-existence. It is not uncommon for a mariner to 
report that, being in the vicinity of a charted rock or shoal, he laid his course 
so as to pass over the geographical position assigned to it with one hundred 
fathoms of line out or with lookouts posted aloft, but was unable to detect any 
evidence of its existence, and that he does not believe, therefore, that the rock 
or shoal exists. . 

It seems necessary, therefore, to inquire into the degree of confidence 
that can be placed in such a piece of evidence of the non-existence of a dan- 
ger, and to establish what probability there would be of finding it under these 
conditions. 

Suppose that A discovers, in the open ocean, a shoal 7 miles in radius, 
and determines the geographical position of its center subject to extreme 
errors of 7 miles in longitude and 7 miles in latitude ; and that B, who is able 
to establish his geographical position within the same limits of extreme error 
as A, attempts to find the shoal again by proceeding to the geographical posi- 
tion assigned to it by A. What is the probability that he will find it ? 

If A, after making the discovery, had revisited the shoal a great number 
of times and had deduced the latitude and longitude of the same spot, under 
the same circumstances, at each visit, the latitudes would all differ from the 
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true latitude, and, likewise, the longitude from the true longitude. If we call 
the differences between the true latitude and the deduced latitudes errors of 
latitude, and lay them off, according to their signs, to the right and left of an 
assumed origin, and then, corresponding to each error as an abscissa, erect an 
ordinate of a length proportional to the probability of that error, these ordi- 
nates and abscissas will be the coordinates of the probability curve. And, 
likewise, if the errors in longitude were found and plotted in conjunction with 
their probabilities, a similar curve would be developed. 

In this investigation the probability curve, ordinarily represented by La- 
place’s formula, y == ce’, will be replaced by two equally inclined straight 
lines 1 2 and AZ’ as shown in figure 1. 


' 


0 Cc B 


Fiaure 1. 


This substitution, which has been employed by Hélie in his Traité de 
Balistique Expérimentale and referred to by Wright in his work on the Ad- 
justment of Observations, causes an appreciable but extremely small error 
which has no practical significance when we consider that, from the nature of 
the calculations about to be made, absolute precision is not to be sought. 

The probability of having an error between OC = x and # + Jz (figure 
1) to the right of the axis OS is equal to sJz. As, in this case, OB and OB’ 
measure the extreme errors, all possible errors are comprised between zero 
and OF, and zero and OL’; and the sum of all the elements which are singly 
represented by s. dz, or the area of the triangle ALL’, should be equal to 
unity which is the measure of certainty. The equation to the straight line A B 
will be, calling # the extreme error 0 and + the intercept on the axis of S, 


But, since the area A BB’ = b X m= 1 or b= = this equation becomes 
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And since x can only vary between zero and m, the probability of having 
an error between z and x + Jz will be: 


(1) 


A = 


The causes which produce the grouping of a number of deduced geographical 
positions around the true one are of two kinds ; one tending to place the de- 
duced latitude to the north or south of the true latitude, and the other tending 
to place the deduced longitude to the east or west of the true longitude. So 
that a particular deduced geographical position /? will be the result of having 
an error (A in latitude aud an error V7 in longitude. 


0B Xx 
Ficure 2. 


The probability that the geographical position deduced by A, upon his 
discovery of the shoal, occupies a certain position with reference to the true 
geographical position of the shoal is, therefore, easily deduced. Through the 
true geographical position of the shoal let two rectangular axes, O.V and ( Y, 
be passed as shown in figure 2. Upon the former conceive errors in longitude 
to be measured, and upon the latter, errors in latitude. The position /’, of 
which the coordinates are 7 and y, results from the concurrence of two condi- 
tions, the error of 1 miles in longitude and the error of y miles in latitude. 
The probability p, of an error between 7 and v + Je is, as shown by equa- 
tion (1), 


and, in the same manner, the probability p», of an error between y and y + Jy 
will be 


= — dy. (2) 
In these formulas, ™ and 7 represent respectively the extreme errors in longi- 
tude and latitude in miles. 

The probability p of having, at the same time, the error x and the error y, 
or of deducing the geographical position /? as the position of the shoal, will 
be the product p, p,, or 

(m — xr) (n — y) Ay, (3) 


nen 


\| 


166 LITTLEHALES. IMPROBABILITY OF FINDING SHOALS. 


an equation in which 7 can vary from zero to m, and y from zero to x. It is, 
therefore, applicable to the first right angle of the axes OX and 0)’, but, in 
order to make it applicable to other quadrants, it is only necessary to change 
the signs of and y. 

Equation (3) then expresses the probability that A’s determination of the 
geographical position of the shoal is in error by x miles in longitude and y 
miles in latitude. 

If the center of the shoal were really located in the geographical position 
assigned to it by A, and B sbould succeed in coming within 7 miles of it, he 
would find the shoal since its radius is 7 miles. 

We have, therefore, as the second step in the solution of the problem, to 
determine what is the probability that B will come within a circular area, 7 
miles in radius, having its center anywhere within the rectangle described 
about the true position of the shoal with sides equal to the extreme errors to 
which the determinations of latitude and longitude by A and B are subject. 

To find the probability, 7’, of coming within any portion of the rectangle 
of extreme errors inclosed by a curve whose equation is y = f(a), it is suffi- 
cient to integrate the expression (3) between limits depending only upon 
y =f (wz), and we shall have, in the first right angle, 


1 ¢ 
P= J dn {(m —x)(n — y) dy. (4) 


\ 


For a circular area of radius 7, we shall have for the first quadrant, 
Vr? — 2? 


f (m — x) dx fin y) dy ; 
0 0 


and for the whole circle, 
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xe x) de J (n — y)dy, 
0 0 
or 


mn (2 Bn | 


The probability that B will find the shoal depends upon the concurrence 
of two independent conditions whose separate probabilities are represented 
by equations (3) and (5), respectively, and is, therefore, equal to P. p, or 


mnt2 38m 4mn mn? 
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Integrating the two expressions which make up equation (3) between the 
limits 2 and x + Jx and y and y + Jy, respectively, the above expression 
becomes : 


mnti2 38m Bn Amn m 


2m n Qn 


which, for 7 = 1 mile, 7 = 2 miles, y = 2 miles, # — 10 miles, » = 10 miles, 


1 ‘ 
6173° That is, under the con- 
ditions stated, B would stand one chance in 6173 of finding the shoal. 


and Jx and Jy each equal to 1 mile, becomes 


U. S. Hyprocrapuic Orrice, Nov. 12, 1895. 
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NOTE ON THE CONGRUENCE 2" — (2n)!/(n!)*, WHERE 2n + 1 
IS A PRIME. 


By Pror. F. Mortry, Haverford, Pa. 


1. There are two ways of integrating cost! xd, one by a Fourier series, 
where we express cos*"t! # in terms of cosines of multiples of z, the other by 
a “formula of reduction.” Taking the integral from x = 0 to # = z/2, and 
equating the two forms of the integral, we have an algebraic identity, which 
will, on a suitable supposition as to the integer ~, yield a theorem as to prime 
numbers, trivial or otherwise. Doing this, we have, first, 


2°" = cos (2n + 4+ (2n + 1) cos(2Qn — 
4. (22 + 1). 2n 3)a + (Qn +1)2n...(n + 2) 


cos (2 — cos 
1 Qn +1 
Orn 2n+1 n 
and second, from the formula of reduction, 
Ar 
2n(Q2n—2)...2 
\ (Gn + 1) Ga —1)...8’ 


so that the algebraic identity is 
2n(2n — 2)...2 
(2m + 1) (2n —1)...3 


1  %&-+1 = +1)2n yn (2n+1)2n.,.(n+2) 


Let 22 + lbea prime p. Let us use the notation = 0, mod where 
a/bis a fraction, to mean that when the fraction is in its lowest terms the 
numerator has the factor c. Then, from the above identity, 


on 2n(2n — 2).. 
(2n + 1) (2a — 


om  2n(2n—2)...2 
(2n — 1) (2n —#)...1 


or 


—(—)"—0, mod p’, 
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MORLEY. ON THE CONGRUENCE (—)" 
or 
on —(—)" 0, mod 7’, 
or 
0, mod p’, (2) 


the left hand member being of course an integer. 

This result is given in Mathews, Theory of Numbers, p. 318, Ex. 16. 

When x = 1, 2, 3, the left hand member of (2) is respectively 18, 250, 
4116, that is 2. 3°, 2. 5°, 2?.3.7% Thus, when p = 5 and 7, the left hand 
member 0, mod p’*, not merely mod p*. TI have to prove that this is so when 
pisa prime > 3. 

2. It is convenient to prove, first, that when p > 3 
(p — 2) 
Using the notation and results of Chrystal, Algebra, Vol. ii, p. 525, let 

(7 +1)(@+ 2)...(¢e+ ...4 

We have 


-0, mod p. (3) 


or since A, », A,; 0mod p, if p > 3, 
+ +...+ (p 1) 0 mod p. (4) 
Now (writing always 2x + 1 = p), 


for each numerator has, and each denominator has not, the factor 2” +. 1. 
Hence, both sum and difference of the expressions 


1,1 1 1 1 


= 0, mod p; hence, also, each expression 0modp. Hence, again, 
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Hence we can remove the even terms from the statement (4), and there re- 
mains (3), 

3. Instead of expressing cos” « as a sum of cosines, let us take the formula 
which expresses cos px as a power-series in cos z. This is shown in treatises 
on trigonometry to be 


(—)" cos pa = p cos x — 


+ (—)" cos” z. 
Multiply by @” and integrate from 0 to § z; then, using (1), 


1 — V*) (p? — 8°) 2.4 
5...p 
Therefore, 
1 12.4. — 1) psi 1 
i. e. (Art. 2) 0, mod p*. 
Therefore 


— 


\n 
1. — 2) 


_0, mod 


whence, as in Art. 1, 


2" — (—)" ony - 0, mod p’, 


where 2” -+{ 1 is a prime p, greater than 3. 


Haverrorp Nov. 5, 1895. 
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EQUATIONS AND VARIABLES ASSOCIATED WITH THE LINEAR 
DIFFERENTIAL EQUATION. 


By Dr. Geo. F. Merzver, Kingston, Ont., Canada. 


My attention was first drawn to this subject by reading a memoir by A. 
R. Forsyth, published in the Phil. Trans. of the Royal Society, Vol. 179 (1888), 
pp. 377-489, and entitled Invariants, Covariants, and Quotient-derivatives 
associated with Linear Differential Equations. 

Also, reference 1uay be made to Professor Craig’s Treatise, Chapter XIIT; 
or my dissertation, entitled Invariants and Equations associated with the Linear 
Differential Equation. 

As much is gained through the adoption of a convenient and brief form 
of notation, I introduce the following : 


dy 


| — (1 2’), 

Ys | 
Yr Ys5 ey Yes Yn | 
=e 


I.—AssoctaTE EQuaTIons. 


A linear differential equation of the nth order 


A —0 n. n —1) Py +. Pay” +. Py 


has n fundamental solutions ¥,, ¥, ..-, Yn. From any two of these ¥,, yg we 

can form the function (7,7, — ¥. Ys), or in the abbreviated notation (4,7), and 

(n —1) 
21 


taken as the solutions of a linear differential equation of the e@—*) th order 
A, 


from the x y’s such functions, in general independent, which may be 


=< 


z 


war 
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Taking any three of the solutions y, we can form functions (y,%,'y,") or (4/77), 


win oe —?) in number, in general independent, and thus forming the 


solutions of a linear differential equation of the oh ne ae =5) th order 


A, 


n! 
4!(n — 4)! 
) which are the solutions of a linear differential equation 
A 3 — 0. 
Proceeding thus, and combining the 7’s five at a time, we obtain the solu- 
n} 
5! (n — 


Similarly, from the combination of four solutions y we obtain 


functions (4,47 


tions of A, = 0 of order 5yt combining them six at a time, we form 
5)! 


n! 
—"_ ; finally, combining the y’s n — 1 

6!(n —6)!° 

at atime, # functions are obtained solutions of a linear differential equation 


the solutions of A, = 0 of order 


of the order 
A,.=°%. 

These equations 1, = 0, A, = 0, A, = 0,..., A,. = 0 are called, re- 
spectively, the first, second, third, ..., (# — 2)th associate equations of the 
original equation A == 0. 

The last, viz, A,_. — 0 has long been known as the adjoint of A = 0, or 
Lagrange’s adjoint equation. Let ..., denote its solutions, where 


| Yoy Ya-ls Yatls CG Yn 


(— 


. . . . . . . 


(n-2 (n—2) (n—2) (n—2) 
y\" eves Ya—1 ’ Yati Yn 


When we consider the determinant 


, , 
Yi Y2 Y 3 Yn 


—2 (n—2) (n—2) n-2 
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we see that it vanishes for all values of / from 0 to » — 2, and as the first 
n — 1 rows give the solutions of the adjoint A,_. = 0, we have 


The determinant also vanishes when we replace the «’s by their deriva- 
tives u,” ... provided + n — 1. Thus we have the more gen- 
eral relation 


yu, + yu +... yu.” = 0, =0,1,...,2 —1;4=0,1,...,n—1) (2) 


These relations between the solutions of an equation and its adjoint be- 
come more interesting when the equation is self-adjoint, i. e. when A, , = 0 
is the same equation as A = 0. 

In the year 1889 I pointed out to Professor Craig and Professor Forsyth 
the following relation between the associate equations: The /th associate of 
the adjoint equation is the same as the (n — 2 — /)th associate of the original 
equation. This is proven in my dissertation and is included in a theorem due 
to Clebsch, published in the “ Abhandlungen der Kon. Gesellschaft der Wis- 
senschaften zu Géttingen, Band XVII,” Ueber die Fundamentalaufgabe der 
Invarianttheorie. 


Equation. 


When A,_, = 0 is the same equation as A = 0 the uw must be linear 
functions of the y, thus: 


Saude, (6 ...,8) @) 
where the a, are constants. 

The condition necessary and sufficient that A = 0 should be self-adjoint 
is that the invariants with odd suffix vanish. This theorem found in an article 
by Professor Brioschi published 1891 was proven by myself two years earlier 
(see Professor Craig’s Treatise, p. 495). 

If in (1) we let 7 = 0 and then substitute the uv from (5) we obtain 


=. 


-te 


To consider more closely the constants a,;,, let n = 3; then A = 0 becomes 
+ 3Py + Py =9, 


(2 3”) _ 
(2 3) 


where 


3P, = — and 7, 
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and 


CY) = + + WY2 + (Gs + WYs + 
+ + G32) + = 0. 


From (3) we obtain 
Uy = = + + 
Us = (YiYo) = Ay Y, + + - 
From and its first and second derivatives we obtain @,,, and by 
solving the equations 


= = WY! + + 
= — BP, = + + 5 


viz. 
(1 2’ 3”) a, = 2(2’ 3”) (2 3’) —3(2 3'PP, — (2 3’F, 
(1 3”) ayy = (2 (3' 1”) + (2 8”) (3 — (28°) (3 1") —3P, (2 3) (31), 
(1 2’ 3”) a, = (2 3’) (V’ 2”) 4+ (2’ 8”) (1 2’) — (2 38”) (1 2”) — 3P, (2 3’) (1 2’). 
As (1 2’ 3”) is constant, and also @,,, @,., @,, 80 must the right members of these 


equations be constant. If in these equations we interchange the suffixes the 
right members remain unchanged in value ; therefore 


Proceeding in a similar way for x = 4, we obtain 
= 4u,76, , (4 3”) — (4 8”) — 6P,(4 3), 
dy, = 4”) — (2 4”) —6P,(2 4), a, = 2”) — (8 2”) — 6P,(3 2’). 


Thus «,, = 0, as the invariant ¢, is supposed to vanish. 
Had we solved for @,,, ¢,,, and a@,, we should have found that ¢,, = — @,,, 
dy = — 4%, and a4, = —a,,. That these results hold for all values of » will 


be seen by the following considerations :— 

lo p.,, op 

= + Yo + y, = 0, 
4 cYn 


but ¢(y’) by equations (2) = 0. 
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Similarly we see that 


=0, (r=0,...,n—2) (4) 


which equations are similar to (1); and as the determinant formed with the y 


4 
and their derivatives is not zero, “A must be proportional to w,; i. e. 


Then it follows that 


+ (Ay + Ay) + (Ag, + + 

+ (Qin + An )Yn = + + + AinYn) 
+ Y, + + (Az + Ay) 


As all the a,,’s cannot vanish we must have / = 0 or 2. For / = 0, 


Gig Oy = = = 9 
and 
ay = — 03 bh, 8) 
For k = 2, 


Giz . 


In order that the determinant of the a’s may not vanish 7 must be even 
when a, = — a, Thus we can announce the theorems, 

When n is odd, the determinant of the a is symmetric. 

When n is even, the determinant of the a is skew symmetric and the equa- 
tion ¢g(y) vanishes identically. 

We then have to do with a linear complex, included in 


wy” + +... (/ = 1, 2, 2) 
or 


(1 2°) + a; (1 +... + a, + ay (2 3) +... 
+ ay +... 4+ G4, (n—1,n)= 0. (5) 
Thus we see that the first associate of a self-adjoint equation is of order 
$¢n(n—1)—1. 
When we take into consideration the relations existing between the asso- 


ciate variables, the application of the possible groups of substitutions shows 
that some of the coefficients a, may vanish. 
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i : By the introduction of a linear substitution we can cause all the a, to 
+e ie 4 vanish except those for which ¢ + 4 =n + 1 and that lie in the diagonal 
’ t , a which is not principal, and these may have the value + 1. For” = 4 (5) 
becomes 
(its) + + = 0. 
Taking A + 6Py’ + + Py =9, where P, = P,, the solu- 
1) tions of the first associate will be 
v, — (12), =(13), »=(14), =(24), =(3 4). 
Mi, '} Usually it is difficult to form the first associate, but the use of the values 
found for simplifies the process. 


(1 2’) 


= 


— 2(1' 2”) —6Pyw + ay, 


ae x 


vill 6 Py’ + = 2”), 
4+ =2 [((1" 2”) + 3 Po" + 6 Pu — a,,)], 


Ses 


— 


A, — — 12 Pw” + + (18P,;" + 36P?2 — 4Pv’ 
+ (6P,;" + 36P,P/ — 2P/)v = 0. 


— 


Forming the invariants 4, and @, for this, we find that they vanish; A, is, 


therefore, self-adjoint. 
It is not difficult to show that A, = 0 has for its first associate an equa- 


tion of which the solations are 


After the substitution we have, 


= (— 1)" Ya”) = » 


a 


where & = n + 1 —a. If in any portion of the plane the y’s have the form 
x”, the exponents belonging to y, being m,, we obtain the following relations : 


i) ] For x = 3, 
i | a = dyy, and m, + m, — 1=m,; whence m, = 1, and m, + m, = 2. 
‘>. ‘ 
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For n = 6, we have 
ye" ys) = ay,, then m, + m, + m, + m+ m,—10=—m,; 


(Ys! Ys" Ya" Yo") = CY» and =m, + m, + m, + m + m —10 = m,; 
similarly, m, + m, + ms, + m;, + m,— 10 = m,; 
m, + m, + m, + m, + m — 10 = m,. 


Adding the first three and subtracting the fourth twice gives 2 (a, + m,) = 10, 
and similarly for the rest ; so that 


m, + m, = m, + m, = m, + m, = 5 


and in like manner 
+ <= mM, + = mM, + =... mM, + =n —1. 


Where y, be of the form e”*, the same relations exist between the m, but 
90 (#2 = 1, 3, 3,4,...). 

In an interesting and exhaustive article in Crelle, Bd. 113, G. Wallenberg 
treats the class of equations known as Fuchsian equations. They have rational 
coefficients, the integrals behave regularly in the neighborhood of the singular 
points, and the roots of the indicial (determinirende, fundamental) equation 
are al] rational numbers. Halphen has shown that an equation may be trans- 
formed into another with constant coefficients if the absolute invariants are 
constant. 

Then for a self-adjoint equation having its absolute invariants constant 
and belonging to Fuchs’ class. Mr. Wallenberg finds the y in the form 


Then there exists » — 2 homogeneous relations of the second order between 
the y ; viz.: 
n = = YsYn-2 = YeYnti-k » 


(k =1, 2, 3,...) 
Yn J 


In this (7) is a rational function which becomes unity so that 


1 


y= =e”, 
Ynt+i-k 


= 
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Then we have 


Substituting the values of the y and taking n = 6, 


1141 
n 
| re 
For xn = 7, 
| 
1 
ry p2 
1 3 


with similar expressions for the other adjoint variables. It is not difficult to 
find them for the equation of the uth order. They serve to verify the preced- 
ing theory. 

When all the invariants vanish we have the case in which A = 0 may be 
reduced to the form y" == 0, so that the solutions are 


x; (4 = 1,%,3,..., 2) 
and we have the parabolic relations 


Yu = Yer Yet - (A = 1, 2,..., — 2) 
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THE CALCULUS OF VARIATIONS. 


By Dr. Harris Hancock, Chicago, Ill. 
INTRODUCTION AND GENERAL OUTLINE. 


Since the time of the Bernouillis, mathematicians have in a greater or 
less degree considered problems which could be solved by methods of varia- 
tions. Euler and Lagrange gave these methods more systematic and compre- 
hensive forms, and founded the calculus of variations on a more scientific 
basis. 

By extending these principles other mathematicians have augmented the 
subject in a wonderful manner ; without, however, avoiding many difficulties 
which arise from want of rigor in the proofs, and from a misinterpretation of 
some of the fundamental conceptions. All these difficulties were removed 
when Prof. Karl Weierstrass, in 1879-80, founded the entire calculus of varia- 
tions on a new basis, free from any objection, and which at the same time is 
more comprehensive in its embrace. 

The writer does not think it out of place to bring these investigations 
before the readers of the ANNALS OF MATHEMATICS, since so little is known of 
this new treatment of the calculus of variations, especially to American 
students; and so through the courtesy of the editor he will give separate 
papers which are, so far as possible, complete in themselves, and which at the 
same time are intended to include the essential parts of the new theory of 
variations. 

These papers are in a great measure abstracts of lectures that were given 
at Berlin by Prof. Weierstrass ; much is also due to Prof. H. A. Schwarz, 
whose lectures on the calculus of variations the writer had the pleasure of 
hearing at Berlin during the summer semester of 1891. 

1. Ln the differential calculus a definite function is given, and a special 
value of the variable or of the variables (if there are more than one variable) 
is sought, for which the function takes the greatest or the least possible value ; 
in the calculus of variations a function is sought, and an expression which 
depends upon this function in a certain known manner is given. A definite 
integral is given, in which the integrand depends upon the unknown function 
in a known manner, and it is asked what form must the unknown function 
have in order that the definite integral may have a maximum or a minimum 
value. 

We treat only real values of the variables. 
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a 2. In order to learn to recognize the general nature of the subject, we shall 

fe first state a few problems which may be solved by the calculus of variations ; 

a so that, while we seek the general characteristies of these problems, we shall of 
o. our own accord come to a more exact statement of the problems which the 

calculus of variations has to solve. 

B 
4 

| 

| 
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Proptem points and with coordinates (a, yo) and (a, 
respectively ave given. Both points lie on the same side of the axis of X in 
the plane wy. It is required to join P, and P, by a curve such that when the 
pline is turned through one complete revolution about the axis of ax, the zone 
generated by this curve may have the smallest possible surface. 

3. To show some of the defects of the old methods we proceed as fol- 
lows: With the assumption that it is possible to draw a curve through the 
two points which satisfies the conditions of the problem, we suppose that two 
points 7” (2, and P” are taken on the curve, and we find another 
point 7”, y) on the curve such that 


” 


We suppose that /?and 7”, Pand /?” are joined together by straight lines ; 
and later we suppose that these three points are taken very close together, so 
that there is a transition from the two straight lines to the curve. The remain- 
ing portions of the curve on the left hand side of /” and on the right hand 
side of /?” are supposed to remain unaltered. 

The portions of surface generated by the straight lines /’?7” and P/??” are, 


respectively, 
| + (de + (y— yf and + + — 
md The sum of these two surfaces of revolution we consider as a function of the 


variable y, and it is required to find when 


4! is a Minimum. 
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In order to have a minimum this expression when differentiated with re- 
gard to y must be zero; i. e. ] 


zy + (y—y (dey? + — yf 
(A) 


y may be determined from this equation, 
y = f (2), say. 
Therefore i 
=f(«# — de), 
and 

= f(@ + Jz). 
Hence, by Taylor’s theorem, 


=f (e — te) = —f (a) de + (de)? 
y = de) =f (x) + (x) de + (2) 4+ .... 


Hence | 
y—y =f (a) de —3 +..., 


—y =F (a) de + +... 


Substituting these values in (A), we have 


+ 

+ 


Dividing through by Jz and making Jr = 0, we have 
1+ — f(x) = 9; (B) 


(dy\*_ 
1+ de, Y de 


Therefore in order to have a minimum value, /(7) or y must satisfy this dif- 
ferential equation; however, when y satisfies this differential equation, we do 
not always have a minimum, as will be shown later. 


tls 
= 


= 


32 


182 HANCOCK. THE CALCULUS OF VARIATIONS. 


Differentiate (B) with regard to 2, and we have 


dy ay 


da’ da dz’ 


or 
dy ad d | 
dx dx da? 
dx’ 
2 
That is, y = ¢ d’y where c’ is the constant of integration. 


dx*’ 


Since y = e*, and y = ¢-° are two solutions of this last differential equa- 
tion, the most general solution is 
x 
y=ce + °, 
where ¢, and ¢, are also constants. 

This last equation is the equation of the catenary. 

4. Thus, by the help of the theory of maxima and minima, we have, it is 
true, come to a certain result ; but, on the other hand, we have yet to ask 
whether this curve gives a true minimum, and owing to the manner in which we 
have come to the result, we have yet to see whether this curve only in a definite 
portion or throughout its whole extent possesses the property required in the 
problem. 

That we are justified in insisting upon this last statement is seen from 
what follows later, where it will be shown that the curve found above satisfies 
the required conditions only between given limits. 

A simple consideration shows that the method we have followed above is 
not at all rigorous ; since it presupposes, which of itself is not admissible, that 
the curve which satisfies the problem is regular in its whole extent, since 
otherwise the portions of curve between the two points (« — Jz, 7’) and (x, y) 
could not be replaced by straight lines joining these two points. 

5. The characteristic difference between problems relative to maxima and 
minima and the problems which have to do with the calculus of variations, 
consists in the fact, that in the first case we have to deal with only a finite 
number of discrete points, while in the calculus of variations the question is 
concerning a continuous series of points. 

If we wish to substitute in the place of the curve first a polygonal line 
and afterwards apply to this line methods similar to those used above, then it 
turns out that, after we have found a line which satisfies all the conditions, it 
is necessary yet to prove that the required limiting transition from polygonal 
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line to curve in reality results in a definite curve which satisfies the conditions 
of the problem. 

The method given above has been chosen to make clear what is in com- 
mon between, as well as the difference between, the theory of maxima and 
minima and the calculus of variations, and we shall now formulate the problem 
in a different manner. 

6. Every limiting transition, as from polygon to curve, is made of itself, 
if we make use of the conception of integration ; since an integral represents 
the limiting value of a sum of quantities which increase following a definite 
law so as to become infinite in number, the quantities themselves becoming 
smaller and smaller in a corresponding manner. 

If we, therefore, define the surface area of the curve y = f(x) which we 
have to find by 


or 


then this integral will have a definite value for every curve that is drawn be- 
tween the points /’, and /), and consequently the problem may be stated as 
follows :— 

Propiem I. y ts to be so determined as a function of x, that the above 
integral shall have the smallest possible value. | 

The solution of this problem will be given later. 

7. As a second problem may be given the problem of the 4rachisto- 
chrone (curve of quickest descent) which may be stated as follows :— 

ProsieM II. Zn a vertical plane a curve is to be drawn from a point A to 
a point B below in such a manner that a material point which is acted upon 
by gravity, and which is compelled to move upon this curve, shall with a given 
initial velocity go from A to Bin the shortest possible time. 

Let the mass of the material point be 1, its initial velocity a, the accelera- 
tion of gravity 29, the time ¢, and the coordinates of A and / respectively 
(0, 0) and (a, 4). 

Let the direction of the positive y-axis be the direction of a falling body 
(gravity), and let the positive z-axis be directed towards the side on which 2B 
lies. Then according to the law of the conservation of energy, 


(dx\? (dy)? __ 


| 

S = 22 yds, 

= J 
2, " 
| 
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or 
(dx |? 
1+ | ay 
+ dy ; 
whence 
, 


dy | 
=f 4gy +a 
0 


Our problem then is to so determine 2 as a function of y that the above 
integral shall have the smallest possible value. 

8. In the two problems given above one of the variables is a one-valued 
function of the other; this is due to the fact that the system of coordinates 
may be so chosen in both cases. Since this is not possible in all cases, it is 
expedient to represent the curve by two equations, that is, to consider # and y 
as one-valued functions of any quantity ¢, where ¢ has only the property, that 
when it goes through all values between two given limits, the corresponding 
point , 7 traverses the curve from the beginning point to the end point, and 
in such a way that to a greater value of ¢ there corresponds a later* point of 
the curve. Hence the integrals of our two problems, which are to have a 
minimum value may be expressed in the form 


t 
Ss = an uV a? + dt, (1) 
at, II 


where 


9. Propiem IIT. Between tivo points on a regular surface f (x, y, 2) = 9, 

a curve is to be drawn so that its length is a minimum, 
Consider the orthogonal coordinates 2, y, 2 of a point of the surface repre- 
sented as one-valued regular functions of two parameters wv, v. If we consider 
these as the rectangular coordinates of a point of the plane, then to every 


P, 


*If¢, < t, and the point 7, corresponds to ¢,, and P, to t,, then a 
P, in reference to ?, is known as a later point, and P, in reference to 
P, as an earlier point. 
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point of the surface there will correspond a definite point of the v-plane, and 
these points in their collectivity fill out a definite part of the plane, which may 
be looked upon as the image of the surface on the plane. To every curve on 
the surface corresponds a curve in this part of the w7-plane and reciprocally. 

Consider, further, v and v as onevalued functions of a quantity ¢; hence 
to every value ¢ there corresponds a point of the wv-plane, and therefore, also, 
in case this point lies in the definite part of the wv-plane, there is a corre- 
sponding definite point of the surface. 

Consequently if 4, and 7, are values of ¢ which correspond to the two fixed 
points on the surface, then the length of any curve which lies between these 
two points is determined through 


| pf ¢ du dv. »(dv\? 
where 
(on)? 
P= | | + + | | 
Lou | Let Lee | 
_ Once , dydy , d2 22 
v= ou on Ou ov ou ov’ 
2 2 { 2 


We have to determine w and v as functions of 7, so that Z is a minimum. 

10. What is common to these three problems is that we have to deter- 
mine two functions of ¢ in such a way that an integral depending upon them 
and of the form 


t, 
=J F(a, y, dt 


to 


may have the smallest possible value. 

Here ¢, and ¢, have fixed values, and /'(, y, 2’, y/) represents a function 
of x, y, x, y’, of which 2’, y’ may be regarded as unlimited variables, while «, y 
are limited to a region which extends, over the whole plane, or over a contin- 
uous part of the same. 

11. The condition that ¢,, 4, should have fixed values is not essential ; 
moreover both end points may move, as in the case of the third problem, if 
we give it the following form: Two curves are given on a surface ; among all 
the possible curves between the points of the one curve and the points of the 
other curve, that curve is to be found which has the shortest length. We are 
accustomed to call this the geodesic distance of the two curves. 
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In order to solve this problem, we must first solve the special problem 
III, since if a curve has the property of being of minimum length such as is 
required above, it must also retain the same property, if we consider the end 
points fixed. Hence, from III the nature of the curve must be determined. 
The variation of the end points give in addition certain special properties, 
which the curve must possess. 

For example, the shortest distance between two curves which lie in the 
same plane is clearly a straight line; through the variation of the end points 
it follows that this straight line must be at the same time perpendicular to both 
curves. 

12. Proptem IV. Essentially different from the three problems given 
above is the following: We are so to construct a closed curve with given 
periphery, that the surface inclosed shall have the greatest possible area. 

Consider and y such functions of ¢, say x(t), y(¢), that for two definite 
values ¢, and 7, of ¢, the corresponding points 2, y of the curve fall together, 
and that, if ¢ goes from a smaller value ¢, to a greater value ¢,, z, y completely 
traverses the curve in the positive direction ; then twice the area of the sur- 
face inclosed by the curve will be expressed by the integral 

ti 
I =f (ay — dt, 


ty 


and the periphery of the curve is determined by means of 


a? + y* dt. 
to 


13. The proposed problem is now as follows : So determine x, y as func- 
tions of t, that the integral I which depends upon them, shall have the greatest 
possible value, while at the same time I, has a given fixed value. 

Problems of this nature are the. most interesting and of the most frequent 
occurrence. They require a treatment essentially different from that of those 
first mentioned. 

These problems are sufficient to give a conception of the nature of the 
problems which are to be solved by the calculus of variations, and with these 
as a basis it will be possible to define the object of the calculus of variations. 

We must yet, however, introduce the fundamental conception of the varia- 
tion of a curve. In former times the calculus of variations was considered 
one of the most difficult branches of analysis ; it was wrongly thought that 
the ground of this difficulty was in the supposed lack of clearness in the fun- 
damental conceptions, especially in the conception of the variation of a curve, 
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while the difficulties that do arise, lie for the most part in quite another diree- 
tion. 

14. In the theory of maxima and minima we say the value of a function 
is, for a definite system of values of the variables, a maximum or a minimum, 
if this value of the function for this system of values is greater or smaller than 
for all the neighboring systems of values. 

We say of a function f (x) of one variable, it has at a definite position 
2 == a,a maximum or a minimum value, if the value for x = a ix respec- 
tively greater or less than it is for all other values of x which are situated in 
the neighborhood \x—a_ < 0 as near as we wish to a. 

The analytical condition that /(2) shall have for the position « = a 


a marimum is expressed by — f(a) < 0; 
a minimum F(x) — f(a) > 0. 
In the same way, we say of a function f(a, ty, ..., #,) of n variables, 


that it has for a definite position x, = dy, Uy... Uy mavimum 
or a minimum, Uf the value of the function for x, = dy, ® = Uy, ..., Ly, = Ay 
is respectively greater or smaller than it is for all other systems of values, 
which are situated in the neighborhood |r. ay) < (4 = 1, 2,..., as 
near as we wish to the first position, 

As here we speak of a neighboring system of values, so also in the caleu- 
lus of variations, we speak of curves which lie in the neighborhood of a given 
curve ; and we require that an integral in the case of a minimum should be 
less and in the case of a maximum greater when taken over the given curve 
than for any of the neighboring curves. 

15. In order to fix the conception of a neighboring curve, and to make 
clear the analogy of the same with the conception of a neighboring system of 
values, let us consider first instead of the given curve a broken line A,A,A, 
... A,, and let us cause the same to slide just a little from its original position. 

Then in the new position every corner /, will correspond to a definite 
corner A, in the old position, and moreover the new position 2, 2,B,... B, 
will be as little different from the old position A,A,A,... A, as we wish, if 
we stipulate that the distance between any two corresponding points A, and /, 
shall be smaller than any quantity 0 where ¢ is as small as we choose. Now, 
by increasing the number of sides, let the broken line pass into the given 
curve, then the points B,B,... 4, will also form a curve which is little dif- 
ferent from the first curve, and which we consequently call neighboring to the 


first curve. 
We can, therefore, say a curve is neighboring to another curve, or exists 


£15 
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out of another curve through a variation as small as we choose, if to every 
point of the latter curve there corresponds a definite point on the former curve 
and also the distance between any two corresponding points is smaller than 0, 
where dis as small as we choose. 

This geometrical conception of a neighboring curve offers no obscurity. 
In a similar manner it is easy to see that for every change of the curve, there 
is a corresponding change of the integral 


(Fo, yw, 


and that this change will be infinitely small, when the second curve is neigh- 
boring to the first. 

This change of the value of the integral must of course be a continuous 
negative one, if the integral is to be a maximum, and a continuous positive 
one, if the integral is to be a minimum. 

16. In accordance with this we may formulate the problem of the caleu- 
lus of variations as follows : 

The quantities 7, y are to be so determined as functions of a quantity 7, 
that when we define a curve by the equations 2 — x(t), y = y(f), and cause 
the curve to vary as little as we choose, the change which in consequence takes 


place in the integral 
t 


ty 
i= F(a, y, x,y) dt 
is continuously positive if a minimum, and continuously negative if a maxi- 
mum is to enter. 

If we consider the problem as proposed in this manner, we have a definite 
problem of the calculus of variations before us, and we have to find strict and 
rigorous methods for the solution of this problem. 

17. Variation of curves. 
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Let y = f(z) be the curve which gives a minimum surface area when we 
rotate this curve about the axis of 2. 

Let 7 be the distance between this curve and any neighboring curve 
mneneal on the y-ordinate, and suppose that 7 is a continuous function of x 
subject to the conditions, that for # = 2, 4 = nollie for 7 = x, 4 = 0; and for 
all other points 7 < yp, where » may be as ane as we deme. 


The integral of any neighboring curve corresponding to (1) is 


(y + + ae. (2) 


Hence the total variation caused in (1) when, instead of y — 7 (7), we take a 
neighboring curve, is 


JIS has always a positive sign since the surface in question is a minimum. 
18. Instead of the one neighboring curve, we may consider a whole 
bundle of such curves, if for 7 we substitute <7, where ¢ is independent of « 
and has any value between — 1 and -- 1. After this substitution (3) becomes 


‘ (dy \* 


and, on the other hand, developing JS by — theorem, 
2 
JIS = cos 1-3? oN ON (5) 


There is no constant term in this last development, since when ¢ is made zero 
in (4) the first and second integrals cancel each other. 
ON is the first variation, 


OS is the second variation, ete. 


Instead of taking 7 a very small quantity, we may take ¢ so small that ey 
is as small as we choose. 


| 
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With Lagrange, writing 7 = dy, it is seen that the total change in y is 
ey edy = Ay. 

Remark. The sign of differentiation and the sign of variation may be in- 
terchanged ; for example the 1st derivative of a variation is equal to the 1st 
variation of a derivative, as is seen by writing 


dy, then = (dy! = (ay). (1) 


Again 7 = dy; change y into y + ey, and consequently 7/ into 7 + ¢7' ; whence 


ay ; (2) 
hence, from (1) and (2), 
d _ afdy) 
de 


19. Returning to (4) and writing y/ = = of = - and expanding the 


expression under the sign of integration 


(y + ey)V1+ (y + 
we have 


Hence, equating the coefficients of the 1st power of ¢ in (4) and in (5) we have 


which is a homogeneous function of the 1st degree in 7 and 7 (7 cannot be 
infinitely large, since then the development would not be necessarily con- 
vergent). 

In a similar manner we may find a definite integral for the 2nd variation, 
in which the integrand is an integral homogeneous function of the 2nd degree 
in 7 and 7 ; similarly for the third variation, ete. 
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